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Stable  Model  Reference  Adaptive  Control 


in  the  Presence  of  Bounded  Disturbances 


Gerhard  Kreisselmeier  and  Kumpati  S.  Narendra 


Abstract 


The  adaptive  control  of  a  linear  time- invariant  plant 
in  the  presence  of  bounded  disturbances  is  considered.  In 
addition  to  the  usual  assumptions  made  regarding  the  plant 
transfer  function,  it  is  also  assumed  that  the  high  frequency 
gain  kp  of  the  plant  and  an  upper  bound  on  the  magnitude  of 
the  controller  parameters,  are  known.  Under  these  conditions 
the  adaptive  controller  suggested  assures  the  boundedness  of 
all  signals  in  the  overall  system. 
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Dynaraik  der  Flugsysteme  D-8031  Oberpfaf fenhofen,  West  Germany, 
January  through  June  1981  and  is  visiting  the  Department  of 
Engineering  and  Applied  Science,  Yale  University. 


1.  Introduction;  A  major  step  in  the  development  of  adaptive  systems  theory  was  the 


establishment  in  recent  years  of  the  global  stability  of  several  equivalent  adaptive 

* 

schemes  [1-3].  When  the  sign  of  the  high  frequency  gain  k^,  the  relative  degree  n 
and  the  order  n  of  the  transfer  function  of  a  plant  with  zeros  in  the  left  half  plane 
are  given  it  was  shown  that  these  schemes  can  be  used  to  adaptively  control  the  given 
plant  in  a  globally  stable  fashion.  The  latter  implies  that  the  parameters  and 
signals  of  the  plant  and  the  controller  are  bounded  while  the  error  between  the 
plant  output  and  the  output  of  a  reference  model  tends  to  zero.  The  question 
naturally  arises  as  to  how  well  these  schemes  perform  in  the  presence  of  external 
disturbances.  In  a  recent  report  [4]  it  was  shown  that  whin  bounded  output  dis¬ 
turbances  are  present  the  parameter  error  vector  <5>(t)  can  grow  without  bound  even 
though  the  state  error  vector  between  plant  and  model  is  bounded.  This  can  be 
attributed  to  the  fact  that  the  overall  nonlinear  adaptive  system  is  only  uniformly 
stable  but  not  uniformly  asymptotically  stable.  Hence  it  became  evident  that  modi¬ 
fications  in  the  basic  adaptive  schemes  are  necessary  when  bounded  external  dis¬ 
turbances  are  present. 

One  such  modification  suggested  in  [4]  resulted  from  the  observation  that  the 
adjustment  of  the  parameter  error  vector  (given  by  4>(t))  is  known  to  be  in  the 
"right"  direction  only  when  the  output  error  e(t)  is  large.  Hence,  it  was  argued, 
if  a  bound  on  the  disturbance  is  known  and  $  is  adjusted  only  when  the  output  error 
exceeds  a  computed  threshold,  the  signals  and  the  parameters  of  the  system  would  be 
bounded.  Such  an  adjustment  corresponds  to  an  adaptive  law  with  a  deadzone  and  in 
[4]  it  is  shown  that  the  above  conclusions  are  indeed  true.  The  presence  of  the 
deadzone  however  results  in  finite  parameter  and  output  errors  even  when  no  external 
disturbance  is  present. 

In  this  paper  an  alternate  approach  is  taken  which  retains  the  potential  of 
obtaining  zero  output  error  e(t)  (and  zero  parameter  error  4> ( t )  if  the  input  is  rich) 
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in  the  limit  when  no  external  disturbance  is  present.  Assuming  that  the  desired 
constant  controller  parameter  vector  9  has  a  norm  less  than  a  known  upper  bound, 
the  adaptive  law  is  suitably  modified  when  Il9(t)||  exceeds  this  bound.  This,  in 
turn,  assures  at  the  outset  the  boundedness  of  the  parameter  error  vector  <p(t) 
despite  the  presence  of  the  external  disturbance.  Demonstrating  that  all  the  other 
signals  within  the  adaptive  control  loop  also  remain  bounded  is  the  principal 
result  of  the  paper. 

Surprisingly  enough,  the  difficulty  in  the  proof  of  stability  does  not  arise 
when  some  of  the  state  variables  grow  rapidly  (e.g.  exponentially)  but  rather  when 
they  grow  slowly  with  time.  In  addition  l!x(t)ll  (where  x  is  the  state  of  the  overall 
system)  need  not  grow  monotonically ,  so  that  the  usual  limiting  arguments  cannot  be 
used  directly.  If  ilx(t)il  is  assumed  to  grow  without  bound,  the  existence  or  an 
arbitrarily  large  interval  [t^t^]  over  which  llx(t)||  is  large,  can  be  established. 

It  is  shown  that  in  such  a  case,  the  effect  of  the  disturbance  is  relatively  small 
and  hence  arguments  similar  to  those  in  the  disturbance  free  case  apply.  This,  in 
turn,  leads  to  a  proof  by  contradiction. 

2.  Structure  of  the  Adaptive  Control  System: 

Let  the  plant  to  be  controlled  be  represented  by  the  equations 


x  =  Ax  +  bu  +  dv, 
P  P  P  P  P  Pi 


y  =  c  x  +  v 
P  P  P  2 


(1) 


where  x  ,u  and  y  are  the  state,  input  and  output  respectively  and  v,  and  v 
P  P  P  1  1 

are  plant  and  output  disturbances.  The  transfer  function  of  the  plant  is  given  by 

T  _i  N  (s) 

C  (Sl-A  )  b  =  k  -S-r-r  =  W  (S). 

P  P  P  P  D  (s)  P 

The  following  assumptions  are  made  regarding  W  (s)  and  the  disturbances  v 

P  1 

and  v2: 


Throughout  this  report  while  representing  a  function  of  time  the  argument 
'  t'  is  omitted  for  the  sake  of  simplicity  of  notation  when  no  confusion  can  arise. 
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(i)  N  (s) ,D  (s)  are  monic  polynomials  of  degrees  m  and  n  respectively, 

P  P 

and  dim(x  )  =  n. 

P 

(ii)  m  and  n  are  known  but  the  coefficients  of  N  (s)  and  D  (s)  are  unknown. 

P  P 

(iii)  N  (s)  is  a  strictly  stable  polynomial. 

P 

(iv)  k  is  known. 


P 

(v)  v^»  are  piecewise  continuous  and  uniformly  bounded  time  functions 

defined  for  all  t  £  R+. 

A  reference  model  is  defined  by  the  equations 


x  =  A  x  +  b  r 

ra  mm  m 

T 

y  =  c  x 

m  mm 


(2) 


and  its  transfer  function  is  defined  as 


T  -1  1 

c  (sI-A  )  b  =  k  — f- 
m  m  m  m  D  (s) 
m 


A  k 

=  f  W  (s) 

k  m 
P 


It  is  assumed  that: 

(i)  D  (s)  is  a  monic,  strictly  stable  polynomial, 
m 

*  A  * 

(ii)  the  degree  of  D  (s)  is  n  =  n-m,  and  dim(x  )  =  n 

m  m 

and  (iii)  r  is  a  piecewise  continuous,  uniformly  bounded  reference  signal. 

The  objective  is  to  control  the  plant  in  such  a  fashion  that  the  output  error 

between  the  plant  and  the  model,  i.e.  e,  =  y  -  y  ,  as  well  as  all  the  state 

1pm 

variables  remain  uniformly  bounded.  This  objective  applies  both  to  the  case  of 

model  following  (r  ?  01  as  well  as  state  regulation  (r  ;  0,y  =  0). 

m 

To  meet  the  control  objective,  a  controller  described  by  the  equations 


v,  =  Fv  +  gu  (3a) 

1  1  p 

v  rt  =  Fv  +  gy  (3b) 

2  2  p 

T 

u  =  v  9  +  —  r  (4) 

P  kn 


is  set  up,  where 
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!  i  .  I  I,  , 
v  =  [V1»V2J  and 

(i)  dim(v^)  =  dim(v2)  =  n. 

(ii)  F  is  an  arbitrary,  strictly  stable  matrix 

(iii)  (F,g)  is  a  completely  controllable  pair 
and  6  is  a  vector  of  controller  parameters  to  be  adapted. 
The  equations  for  the  adaptive  scheme  are: 


r, 

S 

01 

01 


1 

2 

1 

2 


A  C  +  b  u 

ml  m  p 

A  t.  +  b  y 

m2  m  p 


Foi  +  gc  C 
1  ml 


FU>2  +  gCmC2 


(5a) 

(5b) 

(6a) 

(6b) 


0  =  -r 


T  T 
oifoi  0  -  c  £  +  y  ] 


m 


T 

1  +  x  x 


^ —  ref(e) 


f  (©)  = 


(1  -  llell/lle  ii  )  if  II  e II  >  lie  II 

max  max 

0  elsewhere 


T  A  .  T  T  T-i  rT  A  .  T  T,  T  A  .  T  T,  , 
where  x  =  [v  ,£  ,ui  J  ,  £  ^l’^2  ’  W  "*  Wl’U2  and 


(8) 


(9) 


(i)  T  =  T  >  0  is  an  arbitrary  gain  matrix 

(ii)  ||e  ||  is  a  known  upper  bound  on  the  norm  of  the  (unknown) 
max 

* 

matching  controller  parameter  vector  0  . 

Since  v  is  equal  to  the  state  of  a  nonminimal  representation  of  the  plant  (except 


for  exponentially  decaying  initial  conditions  and  the  bounded  effect  of  v^)[l], 

the  controller  structure  allows  the  generation  of  arbitrary  state  feedback. 

•k  k 

Therefore  a  parameter  vector  0  exists  such  that  if  0(t)  =  0  then  the  poles  of 
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the  closed  loop  system  are  the  eigenvalues  of  K  and  the  zeros  of  D  (s)N  (s). 

m  p 

Therefore  the  transfer  function  between  the  reference  input  r  and  the  output 

of  the  plant  y  becomes  k  W  (s) ,  as  desired, 
p  mm 

The  block  diagram  of  the  control  loop  together  with  the  disturbances 

* 

and  v^  is  shown  in  Figure  la.  If  0  =  8  +  $,  where  <$  denotes  the  parameter  error 

vector.  Figure  la  can  be  transformed  as  shown  in  Figures  lb-le  to  yield 
an  error  model  with  a  bounded  disturbance  v  at  the  output.  From  Figures 
lf-lh  it  follows  that 

w(t)T0  -  c\(t)  +  y  (t)  -  v(t)  =  0  ,  (10) 

m  1  P 

v  being  a  bounded  disturbance  due  to  an<^  condition  transients. 

T  A  T  T 

Defining  £  =  [x  ,x  ]  and  using  (10)  in  (8)  the  overall  adaptive  control 

P 

system  may  be  rewritten  in  the  form 

k  =  AC  +  b<j>T(t)v(t)  +  u(t)  (11) 

$  =  -r  it li-t-.y-lt)  1  _  ref(e)  (12) 

1  +  X  (t)x(t) 

where  A  is  a  strictly  stable  matrix  and  p(t)  is  a  uniformly  bounded  input  vector  due  to 
and  r.  The  objective  then  is  to  show  the  global  stability  of  the  overall 
system  (11)— (12) . 

3.  Preliminary  Analysis: 
a)  Adaptive  Law: 

Let  V  =  j  <t>Tr  Evaluating  its  time  derivative  along  the  trajectory  of  (12) 


gives 


T  T 
fr  ui[ui  $  + 

T 

1  +  x  x 


-  4>T9f  (©) . 


pTe  =  (e-e*)Te  5  eTe  -  lle*li  -  lie  l 


From 


-6- 


T  T  4 

and  the  definition  of  f(0)  it  follows  that  <J>  0f(9)  >  0.  Moreover,  <J>  0 f (6)  'v  |UII 
*  T  1/2 

for  IUI|  >>  ||  0  ||  •  Since  u/(l  +  x  x)  is  uniformly  bounded,  this  implies  that 

max 

•  »  f 

V  is  negative  definite  for  all  V  >  V  and  some  V  <  "•  Therefore  V  is  uniformly 
bounded  and  hence  <J>(t)  is  also  uniformly  bounded.  Together  with  (11)  this 
implies  that  £(t)  cannot  grow  faster  than  an  exponential. 

Taking  the  integral  of  (13)  we  obtain 


T  2 
Um) 

T 

1+X  x 


+  4T0f(0)}dT  <  V(t  )  -  V(t-)  +  f  dt 

1  t  J  1+X  X 


5  \  +  A2  *  7 . -TTI72- 

1  t  J  (1+X  x) 


where  A^  and  A^  are  appropriate  positive  constants. 

In  the  absence  of  external  disturbances  the  right  hand  side  of  (15)  reduces 

T  T  1/2 

to  a  constant.  Setting  t^  =  ®  it  can  then  be  concluded  that  w/ (1+x  x)  -*-0 

T  T  1/2 

as  t  +  00 ,  which  in  turn  gives  4>  v/(l+x  x)  -*  0  and  e  +  0  as  t  [5]. 

In  the  present  situation  with  disturbances,  however,  the  integral  on  the 
right  hand  side  of  (15)  may  tend  to  infinity  as  t ^  -*•  even  though  x  -*■  <*>  as 
t  -*■  00 .  Therefore,  the  conclusions  to  he  drawn  from  (15)  have  to  be  modified 
suitably.  This  is  done  in  section  4. 
b)  Control  Loop: 

Let  W  =  £TP£,  where  P  is  a  positive  definite  matrix  and  satisfies  the 

T 

equation  A  P  +  PA  =  -21.  Evaluating  W  along  the  trajectory  of  (11)  gives 

•  T  T  T 

»=-(?+  £  Pb(<Tv  +  u) 


*  -  ££[1 


<cV/2  <tV'2  «V/2 


T  T  1/2 

Since  Hull  and  |£  Pb|/(£  £)  are  uniformly  bounded  there  exist  positive  constants 


3  <  co  such  that 
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CT(t)C(t)  }  a2  (17a) 

|<J>T(t)v(t)|/AT(t)C(t)  <  2ei  (17b) 

implies  that  the  bracket  term  in  (16)  is  greater  than  1/2,  i.e. 

W(t)  <  -  KW(t)  (18) 

where  <  =  1/2  1  (P).  Hence  W  decreases  exponentially  at  least  with  rate  k 

max 

T 

whenever  (17)  holds.  On  the  other  hand,  if  S  (t)£(t)  >  1  then  the  term  within 
brackets  in  (16)  is  bounded.  Therefore  a  positive  constant  X  <  00  exists  such  that 

W(t)  <  XW(t).  (19) 

T 

Hence,  whenever  £  (t)£(t)  5  1,  W  can  at  most  grow  exponentially  with  rate  X. 

In  what  follows  we  make  use  of  the  fact  that  v(t)  is  the  internal  state  of 
a  nonminimal  representation  of  the  plant,  except  for  the  bounded  effect  of  the 
bounded  disturbances.  Therefore  positive  constants  c-^.c^  exist  such  that 
llx  II  s  (1+c  ) II v ||  +  c„-  Hence 

p  x  Z 

I!  X  |]  <  lie  1/  <  II  x  II  +  II  x  || 

P 

<:  CjMxll  +  c ^  (20) 

In  view  of  inequality  (20)  II  €  II  and  llx  II,  when  they  are  large,  are  equivalent  as  far 
as  the  arguments  in  the  following  section  are  concerned. 

4.  Stability  of  the  Adaptive  System: 

Based  on  the  description  of  the  plant  and  the  controller  and  the  assumptions 
imposed  on  them  as  listed  in  section  2,  we  can  now  state  our  main  result  as  follows 

Theorem: 

The  overall  adaptive  closed  loop  system  described  by  equations  (11)  and  (12) 
is  globally  stable  in  the  sense  that  for  arbitrary  bounded  initial  conditions  and 
bounded  signal  p(t)  (or  equivalently  v^(t),v^(t)  and  r(t)),  all  states  of  the 
adaptive  system  remain  uniformly  bounded. 
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The  basic  idea  of  Che  proof  of  the  theorem  is  as  follows.  The  parameter 

error  <Ht)  has  been  made  uniformly  bounded  at  the  outset  by  the  introduction  of 

the  term  r0f(0)  in  the  adaptive  law  (8)  (see  section  3).  Therefore,  only  x(t), 

(or  equivalently  £(t))  can  grow  without  bound.  If  x(t)  grows  without  bound  then, 

since  it  cannot  grow  faster  than  exp(Xt),  it  assumes  values  greater  than  a 

constant  a  over  an  interval  of  time  of  length  a,  where  a  is  arbitrarily  large. 

T  2  T  T 

From  (15)  it  follows  that  (<J>  w)  /( 1+x  x)  and  q>  0 f  ( 0 )  become  less  than  some  t 

T  T  1  /  2 

and  hence  (Ref.  Appendix)  that  |  <f  vj/(l+x  x)  becomes  less  than  most  of 

the  time  within  this  interval.  The  latter  implies  (equations  17b  and  18)  that 
Ux( t) ||  decreases  at  the  rate  exp(-Kt)  most  of  the  time,  and  can  increase  at  most 
as  exp(Xt)  the  rest  of  the  time  within  the  interval.  Since  the  time  interval  over 
which  x(t)  decreases  is  large  compared  to  the  interval  on  which  it  increases  it 
follows  that  ||x(t)||  will  assume  a  value  less  than  a  on  the  interval,  which  contra¬ 
dicts  the  original  assumption. 

Proof: 

Let  us  assume  that 

lim  sup  II  x(t)II  =  »  .  (21) 

t-*»  T<t 


Then,  since  ll  x  II  can  grow  at  most  exponentially  there  exist  monotonically  increasing 

sequences  {t.}>{a.}  with  lim  t.  =  °°,  lim  a ,  =  00  such  that 
l  l  l  ,i 


ll  x(t^)  II  =  (22a) 

!!x(t)l!  >  a.  for  all  t  e  [t.,t.  +  a.].  (22b) 


Using  (22b)  in  (15)  we  obtain 

Xf  {  I;/?  +  ^(eMdt  ?  A  +  A  =  A 

tfJ  (1+x  x)  1  1 

•  T  1/2 

Since  ( t )  and  the  time  derivative  of  w/(l+x  x) 


,  (i  5  0)  .  (23) 

are  uniformly  bounded  (see 


Appendix),  there  exists  a  constant  C,  such  that 
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T  2 
(4>  u) 

T 

1+x  x 


(24a) 


{ 4>Tof (e) }  *  c 


(24b) 


Let  e  €  (0,C]  be  arbitrary  and  consider  i  >  i  such  that  a  5  2CA/c  .  Then 

i 

3 

the  interval  T  A  ft  ,t  +  a  ]  can  be  expressed  as  the  union  of  N  =  2CA/c  disjoint 
i  —  l  i  i 

sub  intervals ,  each  being  of  length  At^  greater  than  or  equal  to  one. 

If  in  any  of  these  subintervals  in  T  there  exists  an  instant  of  time  such  that 
T  2  T  T 

(<j>  cj)  /(1+x  x)  +  $  0 f ( 6)  :>  e  then  (24)  implies  that  this  subinterval  contributes  to 

2 

the  integral  on  the  left  hand  side  of  (23)  an  amount  of  at  least  e  /2C.  Since  the 

2  2 

integral  is  bounded  by  A  there  can  be  at  most  A/(e  /2C)  =  2CA/e  of  these  sub¬ 
intervals.  Let  the  set  of  all  such  subintervals  be  denoted  bv  T . .  •  The  worst  that 

ll 

can  happen  on  T  is  that  W  increases  according  to  W  =  AW. 

3  2 

Let  T  .  denote  the  set  whose  elements  are  the  remaining  2CA(l/e  -  1/e  )  sub- 
12 

T  2  T  T 

intervals  so  that  T^  =  T-qu^12*  ^ave  ^  /(1+x  x)  +  <p  0f(0)  <  e  when  t  £  T  . 

T  T  1  '  2 

It  is  shown  in  the  Appendix  that  this  implies  |  <p  vj/(l+x  x)  S  h(e)  ,  where  h( • ) 

is  a  continuous  function  with  lim  h(e)  =  0.  Hence,  if  t  is  sufficiently  small, 

then  we  have  W  <  -<W  on  T  . 

'  i2 


Therefore 


A  .At  . 

W(t.+(j+l)At.)  <  W(t ,+jAt  )  •  e  J  \  (0  <  j  S  N-l ) 

1  1  11  »  V  J 


where  A  =  -<  on  T  and  A  =  A  on  T  .  This  impl ies 

.1  il  j  12  1 


W(t.+a.)  <  W(t.)  exp{~ —  [-n'(.'-t)  +  At;)). 
ii  l  c  3 

T  T 

Using  £  =  [x  ,x]  and  (20), (26)  we  finallv  obtain 

P 


x(t.+a.)||  <  ||  £(t.+a.)|| 

ii  ii 


A  (P) 

<  IK  (t . )  H  •  — m-a-—  -  -  ■  exp{~  [-<  +  e(A+<) 

1  l\nin(p)J  c3 
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f  a.[c  + 
1  1 


<y3J 


A 

max 


min 


(P) 

Tpy 


•  exp{~ r  [-<  +  e(A+ic)]}.  (27) 

£  J 

Obviously,  the  right  hand  side  of  (27)  can  be  made  less  than  a.  if  e  is 
sufficiently  small.  But  this  is  a  contradiction  to  (22b).  As  a  consequence,  the 
assumption  (21)  was  wrong,  i.e.  x(t)  is  uniformly  bounded,  and  so  is  £(t)  according 
to  (20).  □ 


5.  Conclusions: 

The  model  reference  adaptive  control  problem  in  the  presence  of  bounded  dis¬ 
turbances  is  considered  in  this  paper.  The  principal  result  of  the  paper  is  that 

&  ic 

if  an  upper  bound  U 0  II  on  the  norm  of  the  unknown  controller  parameter  vector  6 

max 

* 

is  known  and  the  adaptive  law  is  suitably  modified  when  I!  0  II  ^  II  0  II  »  then  in 

max 

spite  of  the  disturbances  all  parameters  and  signals  in  the  adaptive  loop  remain 
bounded.  In  the  absence  of  disturbances  the  modification  is  of  no  significance 
and  the  output  error  e(t)  between  the  plant  and  the  model  goes  to  zero  as  t  goes 
to  infinity.  This  implies  that  the  scheme  suggested  retains  the  properties  of 
earlier  schemes  and  is,  in  addition,  robust  with  respect  to  bounded  disturbances, 
which  makes  it  suitable  for  use  in  practical  applications. 

As  mentioned  in  the  introduction,  an  alternate  resolution  of  the  problem  is 
that  if  an  upper  bound  on  the  magnitude  of  the  disturbance  is  known,  then  a  suitably 
designed  deadzone  in  the  adaptive  law  also  guarantees  the  boundedness  of  all 
parameters  and  signals  in  the  adaptive  system  [4].  In  the  latter  approach  the 
descent  property  of  the  parameter  error  vector  ij>(t)  is  always  retained,  whereas 
the  output  error  e(t)  is  only  assured  to  lie  within  the  deadzone  as  t  -*•  <*>,  even 
in  the  absence  of  disturbances.  From  this  it  is  evident  that  the  two  approaches 
behave  quite  differently  both  in  the  presence  of  disturbances  (in  which  4>(t)  may 
not  go  to  zero  in  the  scheme  suggested  in  this  paper)  and  in  the  absence  of  dis- 


-11- 


turbances.  Therefore,  the  analysis  of  their  relative  behavior  deserves  further 
investigation.  Combining  the  two  approaches  appears  to  be  possible  and  also 
interesting  for  practical  applications. 

In  the  adaptive  law,  which  is  considered  in  this  paper,  the  parameter 
adjustment  is  based  on  the  instantaneously  available  information.  Instead, 
all  the  information  that  has  become  available  up  to  time  t  could  be  used  in 
order  to  improve  the  speed  of  adaptation  properties,  as  suggested  in  [5].  The 
adaptive  law  then  becomes 


8 ( t)  =  -r 


tr  w(t)  [u)T(x)0(t)  -  cTC(t)  +  y  (t)]  ,  * 

- - - m_l - E -  e"qU_T)dt  -  T0f (0) 


1  +  X  (t)x(t) 


(28) 


where  q  >  0  is  an  arbitrary  constant.  The  equation  for  the  parameter  error  then 
is  given  by 


$(t)  = 


[oiT(T)4l(t)  +  V(t)  ] 

1  +  xT(t)x(t) 


e-q(t-T)dt  _  r0f(0) 


(29) 


Since  (12)  and  (29)  are  structurally  the  same,  the  stability  arguments  of  this  paper 

can  also  be  applied  to  show  the  stability  for  the  adaptive  law  (29). 

It  is  assumed  in  this  paper  that  the  high  frequency  gain  of  the  plant  is 

known.  This  results  in  a  relatively  simple  proof  of  stability.  The  basic  idea 

of  modifying  the  adaptive  law  based  on  knowledge  of  II 0  II  can,  of  course,  also 

max 

be  applied  when  k  is  unknown.  A  proof  of  stability  may  then  be  obtained  along 
P 

the  same  lines  as  in  this  paper,  although  the  details  may  be  more  involved. 
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(i)  We  first  verify  that  there  exists  a  positive  constant  C  <  »  such  that 


v  <  C. 


II  i/(l+xTx)1/2||  *  C1  (A3) 

II  (jj^iV(l+xTx)1/^2)|  ^  C1  ,  (0  $  i  ><  n  +1)  (A4) 

~  UV^/d+xV72}  <  C1  ,  (0  S  i  s?  n*)  (A5) 

where  denotes  the  i-th  time  derivative  of  u). 

(Al)  holds  by  assumption.  (A2)  was  shown  in  section  3.  (A3)  follows  from 

(i-1')  * 

(A2)  and  (11).  (A4)  is  true  since  C'  '  is  proportional  to  £  for  i  =  l,...,n 

*  T  t 

and  proportional  to  C  and  u  ,y  for  i  =  n  +1.  From  u  -  v  0  +  r  it  follows  that 

P  P  P 

X  1  /  2  i 

u  /(1+x  x)  is  uniformly  bounded  because  so  are  6  and  r.  Since  v  is  equal  to  the 
P 

state  of  a  nonminimal  representation  of  the  plant  except  for  the  bounded  effect 

T  1/2 

of  the  bounded  disturbances  v.  ,v_  it  follows  that  y  /(1+x  x)  is  also  uniformly 

12  p 

bounded.  This  establishes  (A4) .  (A5)  follows  from  (A4)  and  (A2),(12). 

(ii)  It  is  to  be  shown  that 


x(t)  ||  5  1/e' 


+  $T(t)e(t)f(e(t))  < 

l  +  x  (t)x(t) 


'  e  "vh1 


and  ~  ^  1  implies 

|  <t>T(t)v(t)|/[l  +  xT(t)x(t)]172  <  h(e)  t  e  [t  ,t2] 

+  + 

where  { h :  R  -*■  R  }  is  a  continuous  function  and  ouch  that  lim  h(e)  =  0. 
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Since  <j>  0 f ( 9)  >  0,  each  terra  of  the  sum  in  (A7)  is  less  than  e.  Making  use 

of  the  fact  that  f(e)  =  0  for  II 0 II  <  II  0  II  >  we  obtain 

max 


f(e)  <  "t“  =  “t  £ 


T  * 

0  0  0  -  0  0 


0*11  '  II  0 

max 


(A9) 


(12)  together  with  (A7) , (A9)  gives  rise  to 

(  it  it  i 


<  r 


_ _  .  -  H  +  hi 

777  T  ,1/2  ~  T  1/2 

(1+x  x)  (1+x  x) 


+  Hell  •  f (e) 


<  Hr  II  {  JT  +  e3C  +  e/(l|0*ll  -  110*11)} 

i  max 


$  C2  ,  (0  <  e  6  1)  . 


(A10) 


Now  let 


|<j.Ta)(i)/(l+xTx)1/2|  *  e. 


t  e  [t1,t2] 


(All) 


hold  for  some  e  .  >  0  and  0  <  i  <  n  .  This  is  true  for  i  =  0  and  e_  =  /e-  by  the 

l  0 

hypothesis  (A7) .  For  i  +  1  we  obtain  the  following  identity  by  partial  integration 


t+Atr  J  (i+1) 

CO 


,u  T  .1/2 
(1+x  x) 


dt 


_  if  01 


T  (i) 


T  1/2 

(1+x  x) 


t+At  t+At 
t  t 


;T  (i> 

^ —  dt  + 


t+Atr  , T  (i) 
®  01 


T  1/2 
(1+x  x) 


.  T  ”  1/2 
tJ  (1+x  x) 


T. 

X  X 


T 

1  +  x  x 


dr 


(A12) 


where  [t,t  +  At]  <=  [t  .t^].  Defining  C  =  max{l,C  ,C  }  it  follows  that 


t+Atr  aT  (i+1) 

<t>  01 


e  d+A)1/2 


dt 


<  2e  +  At  •  (C  /e~  +  e )  . 


(A13) 
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If  the  integrand  on  the  left  hand  side  of  (A13)  is  equal  to  ei+1  at  time 

*  I 

t  e  [t^,t2],  then  by  choosing  At  =  e  /C  and  t  such  that  t  e  [t,t+At],  we  obtain 

2 

the  result  that  the  integral  in  (A13)  is  greater  than  or  equal  to  At*e:^+^/2  =  ei+i/2C. 
Inequality  (A13)  then  can  be  rewritten  as 


2 

-  2ei+1(c  +  -  *£  C  5  0 


(A14) 


Solving  (A14)  for  the  maximum  possible  we  get 


e  =  (C  7z~  +  e.C)  +  { (C  7z~  +  z  C)  +  4e  CT 
i+1  l  i  i 


(A15) 


In  conclusion,  (All)  holds  for  z  =  /z~  and  e . .  being  defined  recursively 

0  1  x 

by  (A15) .  . 

n  (i) 

Since  v  =  I  d.u  from  (3)  and  (5), (6),  it  follows  that 
i=0  1 


T  1/2 

(1+x  x) 


n  J,  CD 

Z  di  T  ' 
i=0  1  (1+x  x) 


*  I  d  e  -  h(e) 
i=0  1 


(A16) 


Hence  lim  h(e)  =  0,  which  was  the  result  to  be  proven. 

e->0 


1 


